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Low-temperature magnetic properties of both classical and quantum dimerized ferromagnetic 
spin chains are studied. It is shown that at low temperatures the classical dimerized model re- 
duces to the classical uniform model with the effective exchange integral Jo = J(l — S 2 ), where 8 
is the dimerization parameter. The partition function and spin correlation function are calculated 
by means of mapping to the continuum limit, which is justified at low temperatures. In the con- 
tinuum limit the calculation of the partition function and spin correlation function is reduced to 
the eigenvalue problem of quantum rotator in gravitational field. Quantum model is studied using 
Dyson-Maleev representation of the spin operators. It is shown that in the long-wavelength limit 
the Hamiltonian of the quantum dimerized chain reduces to that of the uniform ferromagnetic chain 
with the effective exchange integral Jo = J(I — 8 2 ). This fact implies that the known equivalence 
of the low-temperature magnetic properties of classical and quantum ferromagnetic chains remains 
for the dimerized chains. The considered model can be generalized to include the next-neighbor 
antiferromagnetic interaction. 
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The Paierls instability plays an important role in quasi-one-dimensional materials. The Paierls metal-dielectric 
transition originated from the coupling between electrons and phonons occurs, as a rule, in organic solids [l|. Such 
transition can take place in the quantum spin chains coupled to phonons as well (so-called the spin-Paierls transition 
(SP)). In particular, the SP transition has been observed in the antiferromagnetic spin chain CuGeO^ 0. Currently 
there is a growing interest to the quasi one-dimensional magnets with ferromagnetic exchange interactions and 
the possibility of the Paierls instability in them is widely discussed |6[. It is argued that the possible real system 
where this instability takes place is monatomic chains of Co on the Pt surface Q- As it is proposed 0, these chains 
have the effectively ferromagnetic spin-spin interaction and very weak elastic constants. Recently, another mechanism 
£Sj ' of the Paierls instability in spin systems has been proposed [H, H| . It is based on the coupling of the spins with an 
£> . electronic subsystem (spin-orbital mechanism). It is expected that this mechanism is relevant to the transition metal 
oxide YVO3. The properties of this compound is described by the spin-orbital model [l(| ■ The mean- field treatment 
of this model leads to the ID dimerized spin model with the ferromagnetic sign of the interaction |6|. 

Therefore, the study of the dimerized ferromagnetic (FM) chains is important from both theoretical and experimen- 
tal points of view. The spin chains with the spin-phonon interaction are often described in the adiabatic approximation 
which is valid if the phonon energy is smaller than the Paierls gap. The Hamiltonian of this ID spin model has the 
form 

<N ; 

H H spin ~\~ -^elastic (1) 

where 

H spm = J 5^(1 - (-1)™<5)S„ ■ S„ +1 (2) 



O 
^1" 



_ Nk62 

■^elastic — ^ (3) 

where S is the spin operator, J is the exchange integral, S is the dimerization parameter characterizing lattice 
distortion, k is the effective elastic constant. 
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There is essential difference between the antiferromagnetic (AF) and ferromagnetic (FM) chains with respect to 
the coupling to the lattice distortions. For the AF model (J > 0) the ground state energy of model ^ is ~ — NS 4 ^ 3 
flll [ and the gain in this energy exceeds the loss in Elastic and the SP transition takes place. Contrary to the AF 
model the ground state of H sp i n at J < does not depend on 8 and the ground state of model ([I} has the uniform 
lattice, (5 = 0. However, as was shown in Ref.|6| the thermal fluctuations can activate the dimerization. It was shown 
in Ref.|6J that the free energy of the FM chain at T > is —NT 3 ^ 2 5 2 and the dimerized phase is stable at finite 
temperature for small enough elastic constant k. Such situation can occur in the system Co chains on the Pt surface 
where the Co atoms can be easily moved on the surface. 

The dimerized FM chain can not been solved by the Bethe ansatz in contrast with the uniform model with 5 = 0. 
The thermodynamics of this model has been studied in Ref.[l2[ using both numerical TMRG simulations and the 
analytical modified spin- wave theory (l3| . In particular, the phase diagram of the model has been determined and 
the behavior of spin correlation functions and the zero-field susceptibility are studied. 

It is interesting to consider the influence of the external magnetic field on the thermodynamics of the dimerized FM 
chain. The aim of the present paper is to study the low-temperature magnetic properties of model ([2]) independent 
of the dimerization mechanism. Therefore, we will consider the model with the Hamiltonian having a form 

H = - j£(l - (-1)«5)S„ • S n+1 -hJ^SZ (4) 

where h is the dimensionless magnetic field and J > 0. 

Before we study this model it is instructive to note the remarkable fact related to the uniform FM chain, 5 = 0. 
It was claimed in Ref.[3] that the normalized magnetization M = (S z ) /s of this model at T — > is the universal 
function of the scaling variable g = s 3 Jh/T 2 , i.e. the universal function M(g) is valid for any spin value s and all 
dependence on s is captured in the scaling variable g. In other words, this function is valid for both quantum and 
classical ferromagnetic chains. 

According to Ref.[l4[ the form of the function M(g) can be determined by the computation of the magnetization 
of the classical FM chain in so-called scaling limit, when T — > and h — > but the value of g is finite. As was shown 
in Ref. [3, HH the calculation of the partition function of the classical model reduces to an eigenvalue problem of a 
single quantum rotator. It allows to use the efficient method to compute the function M(g). As a result the explicit 
expansions of this function in small and large g were obtained 

We will show that the universal dependence of the magnetization holds for dimerized model ((U) as well. The 
corresponding universal function coincides with that found in Ref. [l4| but the scaling variable g is renormalized by a 
simple way to include the dimerization parameter 6, namely g is replaced by 7 = g(l — S 2 ). 

The paper is organized as follows. In Section II we consider the partition function of the dimerized classical 
ferromagnetic chain in the magnetic field. We show that the field-dependent part of the free energy coincides in the 
scaling limit with that for the uniform model with the renormalized scaling parameter. In Section III it is demonstrated 
that the spin-wave expansion of the dimerized quantum ferromagnetic chain reproduces the large 7 expansion of the 
magnetization of the classical model. In Section IV the summary of the results is given and the generalization to the 
model with the next-nearest neighbor antiferromagnetic interaction is discussed. 

II. CLASSICAL DIMERIZED SPIN CHAIN IN THE SCALING LIMIT 

In this section we show that the low-temperature magnetic properties of the classical dimerized spin model are 
reduced to that of the FM spin chain with the renormalized exchange integral. 

It is convenient to represent the Hamilton function of the considered classical model Q in a form 

H = —J\ Szi-X ■ §2i — Jl 2_j Sli ' &2i+l ~~ J]] h ' Si (5) 

where J\ = J(l + 6), J 2 = J(l — 6), and Si are spin vectors of the fixed length s and the magnetic field is directed 
along the Z axis: h = (0, 0, h). 

We represent spin vectors on odd and even sites as follows: 

S2i- 1 = sHi 

S 2t = sHi + srhi (6) 

where Hi are unit vectors and rhi are vector differences between neighbor spins (we will assume \rhi\ to be small at 
low temperatures). 
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Then, the scalar products of spins on odd and even bonds become: 



Sai-i • §2i — s 2 — — (jSai-i — S21) 



2 -s 2 



S2i ' S[ 



2i ' >~>2i+l 



-{n'i ~ mi) 2 



(7) 



where we denoted 



H\ = n i+ i - Hi 



After simple algebra the Hamilton function can be transformed to the form 

\ 2 



H 



h 2 N 



2 Ji + J2 
3 — —2^ 



sJ 2 n', 



2(Ji + J 2 ) 2 
Then, the partition function of the system reads: 



s(Ji + J 2 ) 



Z = 



AT/2 

J^J dfhidfii exp 



i=l 



s 2 JiJ 2 
2(Ji + J 2 ) 



g{m»,ni} 
T 



(8) 



(9) 



(10) 



We stress that up to here we did not do any assumption and Eq.© and Eq. (fT0|) are the exact expressions. Now 
we assume that for low temperature T <C (J% + J2) all vectors m, are small and directed in plane perpendicular to 
the corresponding vectors Hi. Then, we can integrate over vectors m, in the infinite limits. That gives 



e 2T(J 1 + J 2 ) 



2nT 



s 2 (Ji + J 2 ) 



N/2 



N/2 



(Hi ■ n lH 



2s 



(ii) 



The first factor in Eq. fTTj) gives a constant contribution to the magnetic susceptibility ~ 1/J. As will be shown 
below, in the low-temperature limit the main contribution is given by the integral in Eq. (jlip and it is much higher 
(~ J/T 2 ). Therefore, we neglect the first factor in Eq. (fTT|) . The second factor does not influence on the magnetic 
properties of the system, and will be omitted. Thus, we reduced the partition function of the dimerized chain to that 
of the uniform ferromagnetic chain with the effective exchange integral 



Jo 



2JiJ 2 
J 1 + -h 



J(l-5 2 ) 



(12) 



The partition function of the classical FM chain in the low-temperature limit was calculated in Ref.[14( by taking 
the continuum limit of the model. We will follow this method. Partition function in the continuum approximation 
takes the form: 



Z ex. I D [H(x)] exp 



dx 
2a 



s 2 Joa 2 ( dH 



dx 



T 



2hs 

1 

T 



(13) 



where L — Na and we notice that the distance between neighbor vectors Hi and n^+i is two lattice spaces 2a, so that 
the vector Hi corresponds to the vector field H(x) at the point x = 2ia in the continuum limit. 

It is useful to transform Eq. (fT3"|) to dimensionless variables. We rescale the spatial coordinate x = ys 2 aJo/T and 
obtain 



Z a J D [H(y)\ exp \- dy ~ 



where A = LT/as 2 Jo is the scaled system length and 



7 = 



s 3 J h 

J^2 



= g(l-6 2 



(14) 



(15) 



Here 7 and g are the scaling variables of the dimerized and uniform models. 
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To calculate the partition function we utilize the well-known equivalence of the n-dimensional statistical held theory 
with the (n — l)-dimensional quantum field theory. The transition amplitude (or propagator) of a particle located 
initially at n(0) — n*o, and finally at n(A) = n\ takes the form of a path integral 

(Hx\ e- xk |n ) <x " D[n{y)] exp j jf L(n' , n)dy^j (16) 

where H (n) is the Hamiltonian operator obtained by quantization of the Hamilton function H corresponding to the 
Lagrangian L. 

Then, imposing the periodic boundary conditions n\ = no = ft and integrating over n, we represent the partition 
function (|14[) in a form 

Z oc J dn(n\e~ x " \n) (17) 
The quantum Hamiltonian has the form (l4| : 

H = l - - jn z (18) 

where I is an angular momentum operator. Hamiltonian (|18[) describes the quantum rotator in the field jn z and 
coincides with the Hamiltonian for the uniform model [14fl with g replaced by 7. 
The corresponding Schrodinger equation in the spherical coordinates has the form: 

1 d 2 cote d m 2 \ , 

Yds 2 ~ —do ~ ^Te ~ 7C0S e ) ^ = £nm ^ nm (19) 

where we used an axial symmetry of the model and introduced the azimuthal quantum number m. 

The exponent of the operator H(n) can be represented using the eigenvalues and the eigenfunctions of the 
Schrodinger equation as follows: 

e -XH{u) = ^ |Vw(n)) e- A£ '"» (i, nm (ft)\ (20) 

ran 

Then, the partition function becomes 

Zoc J dn (n\ e- x " \n) = ( 21 ) 

nm 

In the thermodynamic limit A — > 00 only the lowest eigenvalue £00 (with m = 0) gives contribution to the partition 
function, 

Z -> e~ Aeo ° (22) 
The held-dependent part of the free energy per site is determined by the ground state energy of Eq. (|19|) 



The normalized magnetization M = (n z ) is 



M = -^M (24 ) 
07 



The solution of the Schrodinger equation (|19j) has been found analytically for small and large scaling parameter 
in Ref . [l4| . On the other hand, this equation can be solved numerically for all values of 7 and the magnetization 
curve can be found. It is shown on Fig.l. To demonstrate (5-dependence of the magnetization we represent it as a 
function of the scaled magnetic field g rather than 7. As follows from Fig.l the increase of the dimerization leads to 
the decrease of the magnetization for all values of the magnetic field. 
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FIG. 1: Dependence of the normalized magnetization M on the scaled magnetic field g — s 3 Jh/T 2 for several values of the 
dimerization parameter 8. 



A. Spin correlation functions 



The low-temperature magnetization of the dimerized FM chain coincides with that for the uniform model and only 
the scaling parameter 7 is renormalized. But the spin correlation functions of the dimerized and uniform models are 
different generally. Besides, the longitudinal and transverse correlators are different for the non-zero magnetic field and 
we will consider both types of correlators. At first we consider the correlators (S z S z +2r ) and \SfSf +2r ) = (SfSf +2r ) 
for distances of even number of lattice spacing 2r. In this case the correlation functions do not depend on i and are 
defined by the reduced uniform model (fTTT) : 

(W+2r) = * 2 WO (25) 

The correlator (njn") can be expressed by the eigenvalues and the eigenfunctions of the Schrodinger equation. We 
represent the brief derivation of these expressions on the example of the correlator (rtgrt^). The correlator (ngn*) is 
by definition 



(«0»r) = \ I D[n(y)}n z (0)n z (p)exp ( f 



jf L[n(y)]dy\ (26) 



where 

2Tr 



as 2 Jo 



(27) 



We divide the system on two parts [0,p] and [p, A]. Then the correlator is expressed as an integral of product of 
propagators: 

(n z n z r ) =e Xsm J n z (0)dn(0) / n z (p)dn(p)I(0, p)I(p, A) (28) 
The propagators in the regions [0, p) and [p, A] are calculated with the use of Eqs. (fT6|) and (|20[) : 

7(0, p)= / ' D[n(y)}exp( f P L[n(y)]dy) = £ {n p \^ nm {n)) e^"'" (^ nm (n)\n ) (29) 

and, similarly, 

I(p, X)=^2 (Wnm(n)) e -( A ~P) £ - (ip nm (n)\n p ) (30) 
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where i/j nm and e nm are the eigenfunctions and the eigenvalues of the Schrodinger equation ([19]) . 

Then, imposing the periodic boundary conditions ri\ = ri Q and using the definition n z = J dri \ri) n z (ri\ we obtain 



'o'v/ — 

nrnn' m' 



e-^^'-'-^^-^^-^'-^KVn'm'I^IVnm)! 2 (31) 

ran' m' 

In the thermodynamic limit A — > oo only the lowest level in the region [/?, A] survives, so that n f = m! = and 

The expressions for other types of the correlation functions like (ngn^) has the same form as in Eq. (j3"2")) with 
replacing n z by n x : 

(ng<> = £ |<^oo|^IVw)| 2 exp (- T j^~g° } 2r) (33) 

We expect that the universality in the long distance behavior of the correlation functions holds, so we are interested 
in the asymptotic of the correlation function r » 1. In this limit only the lowest level(s) having non-zero matrix 
element makes contribution to Eqs. (|32[) and (1331) . The operator n z has non-zero expectation value over the ground 
state ("000 1 f^z I "000 ) , which is the normalized magnetization of the system M. Therefore, the main contribution to sum 
(f3"2")l is given by n = to = term and equals M 2 . The decaying correction to this main term is given by the lowest 
excited state with the same azimuthal number to = 0, £io- Thus, the long-distance asymptotic for the correlation 
function (n^n*) is 

(n§<) = M 2 + |W>ook|0io)| 2 ex P (-^-j (34) 

with the correlation length defined by the energy of the lowest excited states as 

£ = g2j ( 1 ~^ 2 ) (35) 
T (eio — £oo) 

The operator n x changes the azimuthal number to, therefore the lowest level for correlator (ngn^) is eqi- This 
implies that the transverse correlation function does not show the long range order and exponentially decays on large 
distances: 

( n o n f) = KV'ooKh/'oi)! exp f - — 

^ = .s 2 J(l-J 2 ) 
T (eoi — £oo) 

The correlation lengths £y and £j_ and the preexponential factors can be found analytically in the limits 7 — > and 
7 — > 00. At 7 — the Schrodinger equation (| 19[) reduces to the equation for the operator of angular momentum with 
well-known spherical eigenfunctions and the spectrum 1(1 + l)/2. So the correlation functions for large r are 

(S?S? +2r ) = (Sf Sf +2r ) = j exp(-2r/0 (37) 



where £11 = £j_ = £ and 



In the limit of high magnetic field (7 3> 1) the lowest eigenfunctions and the corresponding eigenvalues are 

V-oo - 2 7 1 / 4 e^ 2 / 2 , eoo = -7 + V7 
"001 = 7 1/4 0V>oo, £oi = -7 + 2 V7 

0io = O/y^-lMio, £io--7 + 3V7 (39) 
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g 

FIG. 2: Dependencies of | = £,j_T/s 2 J on the scaled magnetic field g = s 3 Jh/T 2 for dimerization parameters 5 = 0, 0.5, 0.8. 



The correlation functions in the limit 7 — > 00 are 



(S?S? +2r ) = s 2 M 2 + — exp(-2r/£u), £|| = - 



5 J(l-5 2 ) 



(SfS? +2r ) = — exp(-2r/£ ± ), U = 



?J(1 - <5 2 ) 



(40) 



According to Eqs. ([57|) and (|4T))) the correlation lengths are changed from £ ~ 1/T at 7 = to £ ~ hr 1 / 2 for 7 — > 00. 
The crossover between two types of the behavior of £ occurs at 7 ~ 1. The dependencies of £_l on g for some values 
of 6 are shown on Fig. 2. 

Now we study the spin correlation function on 'odd' distances. In this case the correlation function (S° l S" +2r+1 ) 
is different for odd and even i. Therefore, we distinguish two types of 'odd' correlators: (S 2 j_ 1 S 2 j +2r ) and 
(S 2 jS 2 j +2r+1 ). The longitudinal correlators (SfS^ +2r+1 ) has non-zero asymptotic M 2 at r — > 00 and the calcu- 
lation of the small corrections caused by the dimerization to this value is not important. On the contrary, the 
transverse correlator (Sf S'f +2r+1 ) decays exponentially and it is more interesting object for the calculation of subtle 
effects like difference on odd and even distances. 

According to Eq.©, the transverse correlator on odd distances is represented as: 



The correlator (n^n x ) was found above in Eq. ([55| . For the correlator (ngm^:) we use the following identity: 



xe 



a{x - y)2 dx = y e- a{x ~ v)2 dx 



Therefore, the integration over m r in the multiple integral 

dmidniriQm x exp 



'0"'r. 



can be transformed as 

J 1 + J 2 



m r exp 



2T 



sJ 2 n' r + h\ 
s(Ji + J 2 ) J 



drh r = 



J 2 n' x 
J1 + J2 



exp 



H{m u fti} 
T 



J 1 + h 
2T 



(41) 



(42) 



(43) 



-»\ 2 ' 

sJ 2 n' r + h \ 
s(Ji + J 2 ) J 



dm r (44) 
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so that 



Using the definition n'* — n* +1 — n% and Eq. ([3l)]) we obtain 

(^•-l^+2r> = ? l^00K|V0l)| 2 e- 2 ^ (l - J2( j~+7 2 /a) ) (46) 

Here we note that for £j_ 3> 1, which is always assumed for low temperature limit, the last factor in Eq. (jlrj]) can be 
expanded and the correlator takes the form 

(^•-l^+2r> = I^OOKIVOI)! 2 e - (2r+1)/ ^ (l - £) (47) 

Similar, for the correlator (Sfj'S'f^r+i) we need merely to exchange J 2 < — > Ji, which gives 

{S^S x 2l+2r+1 )=s 2 \(^o\n x \M\ 2 e-^ +i y^ (l + (48) 

Let us consider an alternation correlation functions [l2j : 

A ± (r) = \(S^ +r )-(S^_ r )\ (49) 
It equals zero for even r. But for large odd r it becomes 



For the small and large 7 we obtain 



A± = s 2 |(VooKl^oi)| 2 ^i e -^ (50) 



The comparison of correlation functions ([37]) and ([51]) for 7 = with those for the quantum dimerized FM model 
obtained in [12| shows that they coincide in the leading terms in T. Therefore, we claim that these correlation 
functions for r 1 obtained for the classical model are valid in scaling limit for the quantum model as well. 

Comparing Eqs. ([47| and (|48|) one can see that the corrections S/£± annihilate each other and gives no contribution 
to the spin structure factor 

which is the sum of these equations. This is valid in the linear in (5/£j_ terms in the low-temperature limit. The terms 
~ (<V£jJ 2 can introduce this dimerization effect into S±(q), but this effect is out of the scope of the used continuum 
approximation. Thus, in the leading term in ST/ J the spin structure factor for the dimerized model coincides with 
that of the FM model with the renormalized exchange coupling. Using Eq. ([33]) for the cor relator (S*S? +r ) we obtain 
the spin structure factor in the form 

S ± ( ) = s4j<yl ~ ^ \^ K^ool^l^" 1 )! 2 ( £ "! Z £ °°) (53) 
T ^ (e„i - e o) 2 + q 2 

with q — qs 2 Jo/T. The dependencies of the normalized spin structure factor S(q) = S±(q)T/ s 4 Jq on q for several 
values of 7 is demonstrated in Fig. 3. 
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FIG. 3: Dependencies of the normalized spin structure factor S(q) = S±(q)T/s 4 'Jo on the scaled wave vector q — qs 2 Jo/T for 
7 = 0,1,3. 



III. SPIN- WAVE EXPANSION OF THE QUANTUM MODEL 



According to the results of Sec. II the magnetization of the classical dimerized model in the scaling limit M(7) 
coincides with that found in Ref. 14 1 . In Ref. [l4[ a method of the computation of this function to an arbitrary 



accuracy was developed. We cite several leading terms of the expansion of M(^) for small and large values of 7 
obtained in Ref. [lij : 



3 7 



M( 7 ) 
M( 7 ) = 1 - 



44 3 

7 + 

135 ' 

1 



752 r 

2835" 
1 



7 < 1 



27V2 1287 3 / 2 512 7 2 



7 > 1 



(54) 



However, it is not clear whether the function M(^f) is universal in the sense that it is valid for both the classical 
and the quantum dimerized FM chains. Below we will produce arguments in favor that such universality is the case. 
At first, we compare the zero- field susceptibility given by 



X(0) 



s 4 J(l-(5 2 ) 



rp2 



dM 
dj 



(55) 



with the asymptotic of x(0) at T — > obtained in Ref.pi 



*<0> = §^ 



for the quantum dimerized FM chain. According to Ref. [l 
S 2 ) 



(56) 



As follows from Eqs. (|54| -([56 |) both expressions for x(0) coincide and function M(j) correctly describes the limit 
7 — > for the magnetization of the quantum model. In connection with Eq. (|55|) we note that it is not applicable for 
the case of full dimerization, 6 = 1. In this case the system consists of decoupled dimers and the susceptibility follows 
the Curie law x(0) ~ 1/T. Thus, our approach is valid when J(l — S 2 ) 3> T. 

Another check of the hypothesis of the universality is the comparison of the spin- wave expansion for the quantum 
model with the expansion of AT (7) for large 7 given by Eq.(|54p. 

The spin-wave expansion is usually carried out by expressing the spin operators using either the Holstein-Primakoff 
or the Dyson-Maleev transformations. Here we use the latter which is 

1 



S: = V2.s(l--a 



2s 



5,, 



2s< 
alar. 



(57) 
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where a+ and a n are the Bose-operators. 

Using Eq. ([57]) we can write Hamiltonian ^ in terms of the Bose-operators. The Bose analog of the spin Hamiltonian 
(j4|) contains terms which are quadratic and quartic in the Bose-operators. The Fourier transform to the momentum 
space operators leads to the Hamiltonian in the form 



where 



H = H + H int - Ns 2 J - Nsh 



Hq = V^[2sJ(l — cosfc) + h]a^a k + 2i8sJ sin ka^ak+K 



(58) 
(59) 



Hint = y^y{k[,k! 1 ,k 2 .,ki)a J [ / a[ / a k2 a kl 5{k[ + k^ - k 2 - k\) 

— fc 1 k 2 

—i5^2 W(k{, k^ 2 ,k 2 , ki)a+ / a^ / a k2 a kl 6(k{ + k^ - k 2 - ki - n) 



v(k{,k£,k2,h) 



W(k{,k / 2 ,k 2 ,k 1 ) 



J r 

4 



cos(fci — k{) + cos(k 2 — k[) + cos(fci — k 2 ) + cos(fc2 — k 2 ) — 2cos(fc() — 2cos(k 2 ) 
'"sin(fci — fc{) + sa\(k 2 — fc{) + sin(/ci — k 2 ) + sa\(k 2 — k 2 ) + 2sin(fc 1 ) + 2sin(fc2) 



(60) 
(61) 
(62) 



Hamiltonian (|59p can be diagonalized by a standard way. Let us rewrite the Hamiltonian Hq using the transfor- 
mation of fc-sums to those over the reduced Brillouin zone and introducing the new Bose-operators a k and f3 k by the 
relation 



a t = u k a++v k p+ 



where |fc| < it/ 2 and 



1 cos k 
2 



iu k B 



(63) 



2e(fc) 

1 cos k 

2 ~ 2e(k) 



e(jfe) = 5 2 )sin 2 k 



Then, the Hamiltonian Hq takes the form: 



where 



H Q = {E a (k)a+a k +E p (k)t3+f3 k } 

|fe|<7r/2 



E a (k) = 2Js[l - e(k)} + h 
E p {k) = 2Js[l + e(k)]+h 



(64) 
(65) 



(66) 



Hamiltonian (|65p describes the non-interacting bosons. At T — > the main contribution to the free energy from 
H is given by the small k region. The expansion for k — > results in 



E a (k) ~ s J(l - 5 2 )k 2 + h 
Ep{k) ~ As J 

so that the thermal occupation numbers of a and /3 particles are 

1 T 



(67) 



n a (k) 



e E a (k)/T _ l 

AsJ 



exp 



E a {k) 
y 



(68) 
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According to Eo.([68|) we can omit in Eq. (|65|) the /? terms giving the exponentially small contribution to the 
thermodynamics at T — > 0. Then the Hamiltonian H takes a form 

H = ^[sJ(\-5 2 )k 2 + h]a+a k (69) 

Eq. (|69|) has a form of the Hamiltonian Ho of the uniform FM model with the renormalized exchange integral 
Jo = J(l - S 2 ). 

Now, let us consider the Hamiltonian Hi nt . First of all, we have to express the operators a k in Eq. ([60|) by the 
operators a k and ft using Eq. ([rJ5|) . As was noted above, for sufficiently low temperatures we can neglect the terms 
in Hi n t containing ft!" and ft operators. Besides, we can replace the Dyson-Maleev vertices V and W by their 
long- wavelength limits. Carrying out some algebra for both terms in Eq. (|60l) we obtain H int in the form 

H int = — - J(l - S 2 ) V] kik 2 a^ / a^ / ak 2 a kl 6(k / 1 + k% - k 2 - ki) (70) 

A remarkable fact is that Eq. ([70)1 is nothing but the quartic in the Bose-operators part of the Hamiltonian of the 
uniform FM chain with the renormalized exchange integral Jo = J(l — <5 2 )- In other words, the Dyson-Maleev vertex 
of the dimerized chain is renormalized one of the uniform model. 

Thus, we established that in the long- wavelength limit, which is justified at low temperatures, the Hamiltonian 
of the quantum dimerized chain ((55)) reduces to that of the uniform FM chain with the effective exchange integral 
Jo = J(l — S 2 ). That is exactly as was found in Sec. II for the classical spin chains. This fact implies that the known 
equivalence of the low-temperature magnetic properties of classical and quantum FM chains remains for the dimerized 
chains. It means that if the universality relatively to the spin value holds for the uniform model then this property 
remains valid for the dimerized FM chain as well. Nevertheless, it is interesting to compare the large 7 expansion for 
the quantum and the classical models. 

The calculation of the spin wave expansion for the free energy and the magnetization in all orders in 1/7 is a 
complicated problem. But the leading terms can be found analytically. In the zeroth order in Hi nt , which corresponds 
to the linear spin-wave approximation, the magnetization in the scaling limit is 

M<°) = - J- (71) 

It is easy to check that the contribution of the first order, vanishes by a symmetry. The two-loop correction 
was calculated in [16| and it is given by 

M {2) = (72) 

I287 3 / 2 v ; 

Thus, the spin-wave expansion of the quantum dimerized FM chain is 



The comparison of Eq.(|54|) with Eq.([73|) shows that both expansions are identical. Though we can not calculate 
the spin- wave expansion in all orders, coincidence of the non-trivial terms in Eq. (|73|) with those for M(p/) is a strong 
argument that the function M("y) gives the low-temperature magnetization of both the classical and the quantum 
dimerized FM chains. 



IV. DISCUSSION 



We studied the low-temperature magnetic properties of the classical and quantum dimerized ferromagnetic spin 
chain. It is shown that at low temperatures the classical dimerized model reduces to the classical uniform model with 
the effective exchange integral Jo = J(l — <5 2 ), where S is the dimerization parameter. The partition function and 
spin correlation function of the classical model are calculated with use of the mapping to the continuum limit, which 
is justified at low temperatures. In the continuum limit the field-dependent thermodynamics depends on one scaling 
parameter 7 = hs 3 J(l — S 2 )/T 2 . The calculation of the partition function and spin correlation function reduces to 
the solution of the Schrodinger equation for the quantum rotator in the "gravitational" field 7. 

We have studied the influence of the dimerization on the magnetic properties of the classical spin model. In 
particular, we have shown that the magnetization decreases with the increase of the dimerization. We found the 
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dependence of the spin correlation functions on both the magnetic field and the dimerization parameter. In contrast 
with the uniform model the correlation functions as a function of the distance r are different for the even and odd r. 
Though the correlation lengths of the spin correlations are the same for even and odd r, the preexponential factors 
are different. 

It was argued in Ref.[14] that the magnetization M of the uniform classical FM chain at T — > is the universal 
function of the scaling variable g — s^Jh/T 2 , i.e. the universal function M(g) is valid for any spin value s and the 
dependence on spin s is captured in the scaling variable g only. It implies that the magnetization curve M(g) is valid 
for both quantum and classical FM chains and can be determined by the computation of the magnetization of the 
classical FM chain in the so-called scaling limit, when T40 and h but the value of g is finite. 

We have shown that this universality holds for the dimerized chain as well. To confirm this fact we studied the 
quantum dimerized spin model with the use of the Dyson-Maleev representation of spin operators. It is shown that 
in the long-wavelength limit, which is justified at low temperatures, the Hamiltonian of the quantum dimerized chain 
reduces to that of the uniform quantum FM chain with the effective exchange integral Jo = J(l — S 2 ). That is exactly 
the same renormalization of the exchange coupling as was found for the classical dimerized spin chains. This fact 
implies that the known equivalence of the low-temperature magnetic properties of classical and quantum FM chains 
remains valid for the dimerized chains. 

The physical reason of the equivalence of the low-temperature magnetic properties of quantum and classical models 
is that the de Broglie wavelength of spin waves As is less than the ferromagnetic correlation length £ [l4j]. Indeed, 
for the spectrum Jk 2 the de Broglie wavelength Xb ~ 1/fc ~ \J J/T, while £ ~ J/T. This implies that the physical 
properties defined by the long-distance asymptotics like magnetization are equal for quantum and classical models. 
But on the short-distances the equivalence failed. For example, the short-distance correlation function behaves as 
(SfSf +r ) ~ s 2 — ar 2 in quantum case [Hj], while the classical model has (SfSf +r ) ~ s 2 — br. 

The considered dimerized ferromagnetic model can be generalized by including in Hamiltonian (j4]) the next-nearest- 
neighbor antiferromagnetic exchange interaction 

J 13 ^S„-S„ +2 , J 13 >0 (74) 

This term leads to the frustration. It is known [l7| that the ground state of the quantum dimerized chain with 
this interaction has the ferromagnetic ground state for a — ^ < (a is the frustration parameter) and the 

1 — H 2 

singlet ground state with the helical spin correlations for a > ° . For the classical model the transition from 

the ferromagnetic to the helical phase occurs at the same value a = 4 . The classical dimerized FM chain with 
the frustration can be studied in full analogy with that for the pure ferromagnetic model. In particular, the low- 
temperature thermodynamics is defined by the solution of the Schrodinger equation (|19[) with the scaling parameter 

7= ^ J( i -*-*,) (75) 

It is easy to check that all presented results are valid for the frustrated model. In particular, the free energy and the 
magnetization are equal to those for the uniform FM chain with the renormalized exchange integral Jo = J(l — 5 2 — 4a). 

Therefore, we believe that the magnetization is described by the universal function, which is valid for both classical 

1— s 2 

and quantum model if the frustration parameter a is not too close to — — . The behavior of the magnetization and 
the susceptibility in the the critical point a = x ~? it is radically different [TH, [Hj]. For example, the zero-field 
susceptibility %(0) ~ T~ 4 / 3 in contrast with T~ 2 behavior for a < 
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